Statistical Mechanics and Dynamics of a 3-Dimensional Glass-Forming System 
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In the context of a classical example of glass-formation in 3-dimensions we exemplify how to 
construct a statistical mechanical theory of the glass transition. At the heart of the approach is a 
simple criterion for verifying a proper choice of up-scaled quasi-species that allow the construction 
of a theory with a finite number of 'states'. Once constructed, the theory identifies a typical scale 
£ that increases rapidly with lowering the temperature and which determines the a-relaxation time 
r a as r a ~ exp(/i£/T) with \x a typical chemical potential. The theory can predict relaxation times 
at temperatures that are inaccessible to numerical simulations. 



Introduction: Among the best studied models of the 
glass transitions are those employing point-particles with 
a soft binary potential. Some repeatedly studied ex- 
amples are the Kobb-Andersen model [l[, the Shintani- 
Tanaka model 0], the Dzugutov model Q and various 
versions of binary mixtures with purely repulsive poten- 
tials, see for example 0, S @|- While easy to simulate 
on the computer, these models are challenging for theo- 
rists due to the fact that it is extremely hard to evaluate 
statistical-mechanical partition-function integrals in con- 
tinuous coordinates. It is therefore very tempting to find 
a reasonable up-scaling (coarse-graining) method that 
would define a discrete statistical-mechanics with parti- 
tion sums rather than integrals, with the sum running on 
a finite number of quasi-species which have well charac- 
terized degeneracies and enthalpies. Indeed, in a number 
of examples in 2-dimensions it was shown that such a 
discrete statistical- mechanics is possible @, H, 0, EH 11 1 
and quite advantageous in providing a successful 

description of the statistics and the dynamics of systems 
undergoing the glass transition. In this Letter we offer 
a general criterion for the selection of up-scaled quasi- 
species and demonstrate it, for the first time, in the con- 
text of a 3-dimensional model system undergoing a glass 
transition. 

Model: We employ here a version of a much studied 
model consisting of a 50:50 mixture of N point-particles 
in 3-dimensions (N = 4096 in our case), interacting via a 
binary potential. We refer to half the particles as 'small' 
and half as 'large'; they interact via a pairwise potential 
U(r t] ): 
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FIG. 1: Color online: Time dependence of the correlation 
functions (f2]) for a range of temperatures (decreasing from 
left to right) as shown in the figure. The inset shows the 
relaxation time r Q in a log-lin plot vs 1/T, compared to an 
Arrhenius temperature dependence. 



The parameter ao is cho- 
0. Below we 
and 
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ra < r c (i,j) 
nj > r c (i,j) 



lates to r c = (a/fS)"-^ t 
sen to guarantee the condition U{r c ) 
use a = 8 and j3 = 6, resulting in r c 
a Q = 0.10546875. 

As in many of these models, one can quantify the slow- 
ing down in the super-cooled regime by measuring a typ- 
ical correlation function. Here we measured the self-part 
of the intermediate scattering function [5j summed over 
the large particles only, 
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(1) 

Here, e is the energy scale and cry = 1.0, 1.2 or 1.4 for 
small-small, small-large or large-large interactions, re- 
spectively. For the sake of numerical speed the potential 
is cut-off smoothly at a distance, denoted as r c , which 
is calculated by solving dU/drij\ rij =r c — which trans- 



F k {t-T) ee ( - ]Texp{zfc • [n(t) - n(0)]} 



(2) 
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In Fig. [T] we show these correlation functions for k — 
5.1(T _1 and for a range of temperatures as indicated in 
the figure. We see the usual rapid slowing down that 
can be measured by introducing the typical time scale 
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FIG. 2: Color online: Temperature dependence of the concen- 
trations of the various quasi-species. Symbols are simulation 
data and the lines are a guide to the eye. 



t q that is determined by noting the time where F^(t = 
r a ;T) = Ffc(0;T)/e = 1/e. The relaxation times are 
shown in the inset of Fig. [T]as a function of 1/T in a log- 
lin plot to stress the non-Arrhenius dependence at lower 
temperatures. 

Statistical Mechanics: Our aim is to provide a sta- 
tistical mechanical theory that captures the structural 
changes upon lowering the temperature such that there 
will pop-up a typical scale that can be used to predict 
the relaxation time r a . To this aim we need to up- 
scale (coarse-grain) from particles to quasi-species that 
can be characterized by their enthalpy and degeneracy. 
Up-scaling can be done in various ways and there is no 
unique algorithm to select a-priori a 'best' up-scaling. 
Here we offer a criterion to validate a chosen up-scaling. 
We choose to work with particles and their nearest neigh- 
bors, where 'neighbors' are defined as all the particles j 
around a chosen central particle i that are within the 
range of interaction r c (i,j). In the interesting range of 
temperatures we find 8 quasi-species with one 'small' cen- 
tral particle and 3, 4 ... 10 neighbors, and 9 quasi-species 
with one 'large' central particle with 6, 7 ... 14 neighbors, 
all in all 17 quasi-species. Other combinations have neg- 
ligible concentration (< 0.5%) throughout the tempera- 
ture range. We denote these quasi-species as C s (n) and 
C(,(n) with s and £ denoting the small or large central 
particle, while n denotes the number of neighbors. We 
measured the mole-fractions (C s (n))(T) and (Ci(n))(T) 
and the results are shown in Fig. [21 
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FIG. 3: Color online: The approximate linear dependence of 
the free energies of the chosen quasi-species on the temper- 
ature. From the slope we read the degeneracy and from the 
intercept the enthalpies (up to normalization), cf. Eq. [4] 
Note that when the free energies are large we do not have 
data: the concentrations become exponentially small and in 
a finite simulation box they disappear completely. 



To decide whether this up-scaling provides a useful sta- 
tistical mechanics we now ask whether there exist free 
energies !F s (n;T) and !Fz{n;T) such that 



(C s (n))(T) 
(Ci(n)){T) 



e -F e (n;T)/T 

2Y 14 e -r e {n;T)/T ■ 
Z— m— 6 



(3) 



The free energies are found by inverting Eqs. ([3]) in terms 
of the measured concentrations. We then plot these 
quantities as a function of the temperature, as demon- 
strated for the present case in Fig. [3J If !F s (n;T) and 
!Ft{n;T) can be well approximated as linear in the tem- 
perature, we can interpret 



T s {n;T) = H s (n) - T In g s (n) , 
Fl(n;T) = H e (n) - T\jig e (n) , 



(4) 



where now the degeneracies g s (n) and ge(n) (read from 
the slopes in Fig. [3]and enthalpies H s (n) and Hi{n) (read 
from the intercepts) are temperature-independent. 
This validates the choice of up-scaling. In other words, 
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the approximate linearity of the inverted free energies in 
the temperature means that we can write the concentra- 
tions as 



(C s (n))(T) 
{Ct{n)){T) 



g s (n)e 



-H.(n)/T 
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f ( n ) e -H.(n)/r ' 



t (n)e~ H ^/ T 



2E:; 6 w(i)e- fl((n)/T 



(5) 



Then we can use these forms also as a prediction for 
temperatures where the simulation time is too short to 
observe the relaxation. The resulting degeneracies g s (n) 
and ge(n) can be easily modeled theoretically, given basi- 
cally by a Gaussian distribution around the most proba- 
ble number n mp of nearest neighbors for small and large 
particles respectively: 

g s (n) » e -[(«-<p) 2 /2^] , <p = 4 .65,ct s 2 = 1.55 , 



(n) « e -[(™-"mp) 2 /2^ 2 ] ) ^ = 7 50j a 2 = 2 Q _ (6) 
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The comparison of the theoretical to the measured de- 
generacies is shown in Fig. [4l upper panel. The same 
figure shows in the middle panel the enthalpies of the 
various quasi-species. One could model the enthalpies 
as a linear function in n. These results are easily in- 
terpreted; we have high enthalpies when there are large 
free volumes (few neighbors). The lowest enthalpies are 
found when there are many neighbors and there is no 
much costly free volume. In other words, at the present 
density and range of temperatures the pV term domi- 
nates the energy in the enthalpy. Using the theoretical 
degeneracies and the measured enthalpies we compute 
the concentrations of all our quasi-species and compare 
them with the measurement in the lowest panel of Fig. 
IH The agreement that we have, especially considering 
the number of quasi-species and the simplicity of the 
theory, is very satisfactory. Notice that the competition 
between degeneracy and enthalpy explains the rather in- 
tricate temperature-dependence of the concentrations of 
the quasi-species, sometimes declining when the temper- 
ature drops, sometime rising, and sometime having non- 
monotonic behavior. 

Prediction of Relaxation Times: Finally, we 
want to connect the structural theory to the dynami- 
cal slowing-down. To this aim we note that there are a 
number of quasi-species whose concentration goes down 
exponentially (or maybe faster) when the temperature 
decreases, and that the relaxation time shoots up at the 
same temperature range. We refer to these quasi-species 
as the 'liquid' ones; in this example the liquid concentra- 
tions are those with small particles with three and four 
neighbors, and large particles with six, seven and eight 
neighbors. We sum up these concentrations and denote 
the sum as (Cii q )(T). The dependence of (Cii q )(T) on 
the temperature is shown in Fig. [5] 
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FIG. 4: Color online: Upper panel: The degeneracies g s {n) 
and gi(n) read from the slopes of Fig. [3] (in circles) and the 
degeneracies according to the gaussian model Eq. (j4|. Mid- 
dle panel: the measured enthalpies. Lower panel: compari- 
son of the measured concentrations of quasi-species to those 
calculated from Eqs. © using the model degeneracies and 
measured enthalpies. Here symbols are data and lines are 
theoretical predictions. 
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FIG. 5: Color online: The temperature dependence of Cii q (T) 
is shown as the upper continuous line. The contributions of 
the various 'liquid' sub-species are shown with symbols which 
are identified in the inset. 
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FIG. 6: Color online: The relaxation time T a (T) in terms of 
the typical scale £(T). We show the excellent fit to Eq. (JSj) 
with /i = 0.04. Note that the intercept at T — ► oo is of the 
order of unity as it must be. 



This concentration is used to define a typical scale 

an 



C(T) ee [pC Uq (T)]-V3 



(7) 



where p is the number density. This length scale has 
the physical interpretation of the average distance be- 
tween the 'liquid' quasi-species. It was argued before 
S US li | that this length scale can be also interpreted as 
the linear size of relaxation events which include 0(£(T)) 
quasi-species. We can therefore estimate the growing free 
energy per relaxation event as AG = n£(T) where \i is 
the typical chemical potential per involved quasi-species. 
This estimate, in turn, determines the relaxation time as 



r a (T) 



(8) 



The quality of this prediction can be gleaned from Fig. [51 
where we can see that the fit is excellent, with fi 0.04. 
The fact that the intercept in Fig. [6J is of the order of 
unity is very reassuring, since this is what we expect when 
T — » oo. 

A few points should be stressed. As we expect (cf. Ref 
in systems with point particles and soft potential, 
there is no reason to fit the relaxation time to a Vogel- 
Fulcher form [l5| which predicts a singularity at finite 
temperature. In our approach we predict that £ — > oo 
only when T — > 0, and there is nothing singular on the 
way, only slower and slower relaxation. At some point the 
simulation time will be too short for the system to relax, 
but we can use Eq. (|8|) to predict what should be the 
simulation time to allow the system to reach equilibrium. 



Conclusions: In summary, we reiterate that point 
particles with soft potential are different from granular 
media or systems of hard spheres which can truly jam and 
lose ergodicity. Point particles with soft potential remain 
ergodic and therefore should be amenable in their 
super-cooled regime to statistical mechanics. To con- 
struct simple, workable statistical mechanics one needs to 
up-scale the system and find a collection of quasi-species 
with well defined enthalpies and degeneracies. In this 
Letter we proposed a simple criterion to validate a choice 
of up-scaling, and demonstrated how, once the structural 
theory is under control, a natural length scale appears 
and can be used to determine the relaxation time, also 
for temperatures that cannot be simulated due to the 
fast growth of necessary relaxation time. The fact that 
the present approach works equally well in two and three 
dimensions provides good reason to believe that it has 
a substantial degree of generality. How to use this kind 
of theory to understand relaxation functions and ther- 
modynamic properties was demonstrated before in two- 
dimensional examples, but there is still a large variety 
of related problem in two and three dimensions that can 
profit from this approach. 
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